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ABSTRACT 

Aims. Deuterated molecules have been detected and studied toward Orion BN/KL in the past decades, mostly with single-dish tele- 
scopes. However, high angular resolution data are critical not only for inteipreting the spatial distribution of the deuteration ratio 
but also for understanding this complex region in terms of cloud evolution involving star-forming activities and stellar feedbacks. 
Therefore, it is important to investigate the deuterated ratio of methanol, one of the most abundant grain-surface species, on a scale of 
a few arcseconds to better understand the physical conditions related to deuteration in Orion BN/KL. 

Methods. Orion BN/KL was extensively observed with the IRAM Plateau de Bure Interferometer from 1999 to 2007 in the 1 to 3 mm 
range. The angular resolution varies from 1'.'8 x U.'& to 3'.'6 x 2'.'3 and the spectral resolution varies from 0.4 to 1.9 km s -1 . Deuterated 
methanol CH 2 DOH and CH 3 OD and CH 3 OH lines were searched for within our 3 mm and 1.3 mm data sets. 

Results. We present here the first high angular resolution (1'.'8 x Qf!8) images of deuterated methanol CH2DOH in Orion BN/KL. Six 
CHiDOH lines were detected around 105.8, 223.5, and 225.9 GHz. In addition, three E-type methanol lines around 101-102 GHz 
were detected and were used to derive the corresponding CH3OH rotational temperatures and column densities toward different re- 
gions across Orion BN/KL. The strongest CH 2 DOH and CH3OH emissions come from the Hot Core southwest region with a velocity 
that is typical of the Compact Ridge (V L sr ~ 8 km s~'). We derive [CH 2 DOH]/[CH 3 OH] abundance ratios of 0.8 - 1.3 x 10~ 3 toward 
three CH 2 DOH emission peaks. A new transition of CH3OD was detected at 226.2 GHz for the first time in the interstellar medium. 
Its distribution is similar to that of CH2DOH. Besides, we find that the [CH 2 DOH]/[CH30D] abundance ratios are lower than unity 
in the central part of BN/KL. Furthermore, the HDO 3i_2 - 2 2 ,i line at 225.9 GHz was detected and its emission distribution shows 
a shift of a few arcseconds with respect to the deuterated methanol emission that likely results from different excitation effects. The 
deuteration ratios derived along Orion BN/KL are not markedly different from one clump to another. However, various processes such 
as slow heating due to ongoing star formation, heating by luminous infrared sources, or heating by shocks could be competing to 
explain some local differences observed for these ratios. 

Key words. Interstellar medium (ISM), ISM: clouds, ISM: molecules, Astrochemistry 



1. Introduction 

Deuterium chemistry in the interstellar medium (ISM) has been 
intensively studied in recent decades. More and more deuter- 
ated molecules hav e been found as well as multiply-deuter ated 
species, e.g., ND3 ( van der Tak et all l2002t iLis et al ., 2002) and 
CD3OH (Pari se et all 2004 ). Deuterium chemistry models have 
also been revised to explain the observations that show a strong 
enhancement of deuteri um-bearing molecule abundances in star- 
forming regions (e.g ., [Roberts et al., 20033 iRoberts & Millaii 
2000: ICharnlev et all 1 1 9971) . compared with the D/H ratio of 
2 - 3 X 10~ 5 in the local interstellar gas (see, e.g.. lLinskv et all 
120061 and references therein). Those strong enhancements (by a 
factor of a few thousand) are seen mostly in methanol, ammonia, 
water, and formaldehyde, leading to the abundance ratios > 0. 1 
compared with their non-deuterated analogs. The formation of 
these molecules can be largely explained by grain-surface reac- 
tions that provide a natural explanation for the abundant doubly 
or multiply deut erated m olecules that ion-molecule chemistry 
failed to predict dTurnerl 1 1990b . 

Many deuterated molecules were first detected toward the 
Orion Becklin-Neugebauer/Kleinmann-Low (BN/KL) region 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany), and IGN (Spain). 



dBecklin & Neugebaueti 119671: iKleinm ann & Lowl fl967l) . one 
of the closest (414 + 7 pc. iMenten et all 120071) and most- 
studied star-forming regions in the sky. For example, deuter- 
ated water HDO and de u terated am monia NH2D were first 
detect ed by iTurner et ail (1 19751) and iRodriguez Kuiper et all 
(fl978l) . and the single-deuterated methanol molecules CH3OD 
and CH2DOH were first detec ted toward the sa me region by 
iMauersberger et al.l (11988b and iJacq et al.l (Il993l) . respectively. 
Although grain-surface chemical models can well explain the 
high abundance of those deuterated molecules in molecular 
clouds, the deuteration branching ratios in the same species 
predicted by models d isagree with observations. For instance, 
ICharnlev et al.l d 1997b showed that the formation of deuter- 
ated methanol on grains based on the addition of H and 
D atoms to CO always leads to [CH 2 DOH]/[CH 3 OD] abun- 
dance ratios of about 3. This prediction is in confli ct with the 
[CH2pOH]/[CH 3 OD] ratio of 1.1 -1.5 observed bv iJaca et al I 
(fl993l) in Orion BN/KL. However. iRodgers & Charnlevl (12002) 
later pointed out that without taking into account other surface 
species such as CO and H2CO and possibly different energy bar- 
riers involved in the reaction scheme, the [CH2DOH]/[CH30D] 
ratio of 3 is an artificial result of the model assumptions. In 
addition, most proposed models that include surface chemistry 
sttongly depend on local environment where temperature and 
gas density play important roles in the chemical reaction rates. 
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Table 1. Observational parameters of the PdBI data sets 



Bandwidth 


Observation date 


Configuration 


Flux conversion 


RMS noise 


^HPBw" 


6\" 


n c 
C'syn 




(GHz) 






(1 Jy beam -1 ) 


(mjy beam 1 ) 


(") 


(km s" 1 ) 


("X") 


PAD 


101.178-101.717 


2003-2006 


BC 


15.8 K 


2.9 


49.7 


1.85 


3.79 x 1.99 


22 


105.655-105.726 


2005 Aug-Nov 


D 


2.9 K 


7.6 


47.7 


0.89 


7.13x5.36 


9 


223.402-223.941 


2003-2007 


BC 


17.3 K 


21.2 


22.5 


0.84 


1.79 x 0.79 


14 


225.805-225.942 


2005 Sep-Nov 


D 


2.9 K 


40.8 


22.3 


0.42 


3.63 x 2.26 


12 


225.990-226.192 


2005 Sep-Nov 


D 


2.9 K 


40.8 


22.3 


0.42 


3.63 x 2.25 


12 



Notes. Primary beam size Channel separation w Synthesized beam size 



This is especially true for the Orion BN/KL region where con- 
tributions from stellar feedbacks (e.g., ultraviolet photons) and 
star formation activities (e.g., outflows/shocks) alter the warm- 
up history of the cloud, involving both grain surface and gas- 
phase chemistry. Hence, it is crucial to investigate the BN/KL 
deuterated methanol distribution with high spatial resolutions so 
that the physical conditions of individual clumps are properly 
constrained and the deuteration ratios can be related to specific 
physical processes. Additionally, deeper insight into the pro- 
cesses at play can be gained by comparing CH2DOH maps with 
another major deuterated species, HDO. 

The massive star-forming region Orion BN/KL is compli- 
cated not only because of the interactions between outflows 
and the ambient material, but also because of its rich and com- 
plex chemistry at the so- called Orion Hot Core and Compact 
Ridge regions (see, e.g., BlakeetaJi [[987). Owing to larger 
single-dish beam sizes in the early spectral line studies of Orion 
BN/KL, molecular line profiles were usually decomposed into 
several components according to their local standard of rest 
(LSR) velocities and line widths. The Orion hot molecular core 
(Hot Core) is usually characterized by its velocity component at 
Vlsr— 5-6 km s _1 and a broad line width of about 5-10 km s _1 . 
It has been identified in the interferometric maps as a strong 
mm/sub mm continuu m emission peak close to the infrared (IR) 
source IRc2 (see, e.g., iGezari eta! .1 fl998l iDownes et all fl98ll 
iRieke et afl , ll973l) . whereas many molecular emission peaks ob- 
served around 8 km s _1 are displaced from the Hot Core by 
about 4" to t he southwes t (Hot Core Southwest, HC-SW; see , 
e.g., works of [F avre et a fl 1201 ll IWang et all l2010t iTang et all 
2010; Friedel & Snyder, 2008, and references therein). 

The 8-9 km s _1 LSR velocity component of the Orion 
Compact Ridge exhibits a relatively narrow line width (3-5 
km s _1 ). The exact location of the Orion Compact Ridge is 
some what ambiguous, but recent interferometric observations 
(e.g., iFavre et all l201lt iFriedel & Snvderl 120081) suggest that 
it is located 10"— 15" to the southwest of the Hot Core (the 
strongest dust continuum peak). It is important to mention that 
the Orion Hot Core is located within the NE-SW d ense ridge of 
the BN/KL region seen in dust continuum emission dFavre et all 
1201 U Fling et alll20 10), and is part of the hierarchical filamen- 
tary struct ure seen on a large r scale in OMC-1 (e.g., mid-7 CO 
images by iPeng et all 120121) . Additionally, the Orion Compact 
Ridge is located at the southern part of this dense ridge, the bot- 
tom part of the V-shaped region see n in many molecul a r lines , 
e.g., CS dChandler & Wood \l99% . SO dWright et all fl 996). 
HCOOCHT ilFavre et all 1201 lh . and NH 3 dGoddi et alll2011h . 

The main goal of this paper is to investigate the deutera- 
tion ratios in Orion BN/KL, and address the possible causes for 
the abundant deuterated methanol in this region. We present the 
first high angular resolution (1'.'8 x 0'.'8) images of CH2DOH to- 
ward the Orion BN/KL region ( 93.1l l. CH3OH maps were also 
obtained from the same data sets with a 3'.'8 x 2'.'0 resolution 



( 9T21 . The CH3OD map and HDO result are shown in 9331 and 
93.41 respectively. We discuss methanol deuteration in 94. II for 
CH 2 DOH and 933] for CH 3 OD. In 9431 our own CH 3 OH data 
are discussed in the light of the spectral line profiles obtained 
at much higher frequencies with Herschel. Comparison of our 
deuterated methanol and methanol maps with our deuterated wa- 
ter maps is presented in 94.41 and water and methanol deutera- 
tion ratios are discussed in 94.51 

2. Observations and spectroscopy 

2.1. I RAM observations 

The data used in this study are part of the large 1-3 mm data 
sets obtained from 1999 to 2007 using toward the Orion 

BN/KL region (see lFavre et alll201 ll for more observational de- 
tails). Four data sets (Table [TJ and five antennas equipped with 
two SIS receivers were used in this study. The quasars 0458-020, 
0528+134, 0605-085, and 0607-157, and the BL Lac source 
0420-014 were observed for phase and amplitude calibration. 
The six units of the correlator allowed us to achieve differ- 
ent bandwidths and spectral resolutions. Using the 30m single- 
dish data (J. Cernicharo, priv. comm.), the missing flux for the 
CH 2 DOH lines around 223 GHz is estimated to be 20%-50%. 
The large uncertainty in this estimate is due to line confusion 
and the difficulty in determining spectral baselines in the 30m 
data. The missing flux in our CH3OH line interferometric ob- 
servations around 101 GHz is estimated to be < 6%. As for the 
HDO, CH 2 DOH, and CH3OD data at 226 GHz, we compared 
our Pd BI HDO line emi ssion at 225896.7 MHz with the 30m 
data of iJacq et all dl990h . and we find that the missing flux for 
this HDO line is about 23% toward IRc2. Since the CH 2 DOH 
and CH3OD emissions are less extended than the HDO emis- 
sion, the missing fluxes of CH2DOH and CH3OD are likely less 
than 23% at 226 GHz. 

The PdBI data were reduced with the GILDAS0 package, 
and the continuum emission was subtracted by selecting line- 
fr ee channels. Our continuum emission images were presented 
in lFavre et al.1 (1201 lh where the H2 column densities of selected 
clumps were also estimated. The data cube was then cle aned 
channel-by-channel with the Clark algorithm dClarkl Il980l) im- 
plemented in the GILD AS package. 

2.2. CH 2 DOH spectroscopy 

CH2DOH is an asymmetric top molecule similar to CH3OH 
where one H atom is replaced by a D atom in the methyl 
group (-CH3). Therefore, the threefold symmetry (C3V) of the 

1 The IRAM Plateau de Bure Interferometer. IRAM is supported by 
INSU/CNRS (France), MPG (Germany) and IGN (Spain). 

2 http://www.iram.fr/IRAMFR/GILDAS/ 
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Table 2. Spectroscopic parameters of the detected lines 



Frequency" 


Transition 






(MHz) 




(K) 


(D 2 ) 


CH 2 DOH 


105806.4100 


Hi, 10 - llo.n e\ 


159.4 


6.11 


223422.2629 


52,4 - 4 2 ,3 e 


48.4 


5.15 


223616.1420 


54,2 _ 4 4 i, 544 - 4 4 e Q 


95.2 


2.15* 


223691.5380 


53,3 - 4 3 ,2 e 


68.2 


3.96 


223697.1880 


53,2 - 4 3 ,i e 


68.2 


3.96 


225878.2540 


3 1,3 - 2o,2 Oo 


35.7 


2.28 


CH3OH 


101185.4530 


6_ 2 - 61 E 


74.7 


0.021 


101293.4150 


7-2-7! E 


91.0 


0.046 


101469.8050 


8_ 2 - 81 E 


109.6 


0.091 


CH3OD 


226185.9300 


5_, - 4_, E 


37.3 


3.4 


HDO 


225896.7200 


3 1,2 _ 22,1 


167.6 


0.69 



Notes. <a) Line frequencies from this work (W Line strength of each 
single line 

methyl group in the internal rotation of CH3OH is broken, 
and the CH3OH symmetry states A, El, and E2 become eo, 
e\, and o\ in CH2DOH, which are more dif ficult to deal with 
theor e tically and experimentally (see, e.g., lLauvergnat et al.L 
2009; iMukhopadhvavetall 120021) . leading to only few pub- 
lished studies of accurate line frequency measurements to date. 
The CH2DOH substate symmetry can be characterized as even 
(e) or odd (o) according to the plane of symmetry defined by 
the DCO or COH plane. For example, the minimum potential 
energy for CH2DOH occurs in the "trans" substate where the 
D atom is opposite to OH, and the secondary minimum oc- 
curs at the "gauche" positions where the D atom is located at 
-120 degrees from the "trans" position (Mukhopad hvav et al.L 
120021) . Therefore, the ground torsional states of CH2DOH con- 
tain three substates, eo ("trans"), e\ (symmetric "gauche"), and 
o\ (antisymmetric " gauche") with an increasing energy, similar 
to eth anol (see, e.g.. lPearson et all [1997: Mukhopad hvav et all 
12001 . The selection rules for the same symmetry states (e.g., 
eo - eo) allow a- and b-type transitions, and c-type transitions 
are allowed for involving different symmetry states (e.g., o\-e\) 
in CH2DOH. The CH2DOH molecular parameters are listed in 
Table [2] The CH2DOH molecular parameters were provided by 
one of us (B. Parise, see Appendix fAl. and the CH3OH molec- 
ular par ameters were ta ken from the JPL database^ (see calcula- 
tions of lXuetalll2008l) . 

3. Results 

3.1. CH 2 DOH 

Six CH2DOH lines have been detected toward Orion BN/KL, 
and three of them are slightly blended with other molecularlines. 
Four CH 2 DOH lines are a-type (AK a = and AK C = 1) R- 
branch (AJ = 1) transitions with one b-type (AK a = 1 and 
AK C = 1) Q-branch (AJ - 0) transition and one b-type R-branch 
transition. The CH2DOH line parameter measurements toward 
Orion BN/KL are summarized in the appendix (Table iBTTT i. 

Emission peaks: The CH2DOH emission is mostly present in 
the southern part of the Orion BN/KL region, or the Compact 



Ridge region with a characteristic Vlsr of » 8 km s . For ex- 
ample, the CFLDOH velocity channel maps in Figures Q] and 
IB. II show one strong peak (deuterated methanol 1, dM-1) and 
two weaker peaks (dM-2 and dM-3) in the 7-9 km s _1 velocity 
range. 

There is no clear CH2DOH detection at the Orion Hot Core 
itself (the strong dust continuum peak ), except for the weak 
emission close to the source / position dMenten & ReidL [1995: 
iGarav et all 1 1987b shown in the CH 2 DOH 5 2 ,4 - 4 2 , 3 e line 
(see Fig.Q] Vlsr=5-7 km s _1 ). Since the feature associated with 
source / is only seen in our strongest CH2DOH line, it is not 
clear if the signal-to-noise ratio (S/N) is too low to detect this 
feature in other transitions, or if this feature is simply caused by 
contamination of other molecular lines from the region close to 
the Hot Core. Nevertheless, there is no indication of other molec- 
ular lines at this frequency. 

The CH2DOH spectra toward the three main peaks are 
shown in the appendix (Figure IB.2t . The average Vlsr mea- 
sured at the peak emission channel is 7.7 + 0.4 km s" 1 for dM- 
1, 7.8 ± 0.4 km s" 1 for dM-2, and 7.8 ± 0.4 km s" 1 for dM-3. 
The FWHM line widths estimated by fitting a Gaussian profile 
are 2.1 + 1.2 km s" 1 and 2.3 + 0.9 km s" 1 for dM-1 and dM-2, 
respectively. The weak dM-3 source exhibits a line width of ap- 
proximately 3.7 + 0.5 km s _1 , which is broader than the dM-1 
and dM-2 line widths; this is likely because of the blending of 
two or more velocity components (see below). 

Additionally, dM- 1 has a counterpart in the methyl formate 
emission (MF-2) observed bv lFavre et all (1201 ll) . and dM-2 and 
dM-3 are nearly coinciding with MF-3 and MF-1, respectively. 
It should be noted that MF-1, t he strongest methy l formate emis- 
sion peak in Orion BN/KL (IFavre et all [201 1), sh ows weak 
CH2D OH emission in our data (dM-3). As shown bv lFavre et all 
d201 ll) . MF-1 exhibits two velocity components at 7.5 km s _1 
and 9.2 km s , and these two components are only seen in the 
225878.3 MHz CH 2 DOH line toward dM-3 where our spectral 
resolution is higher (0.42 km s _1 ). Furthermore, dM-2 and dM-3 
are close to the CH3OH emission peak in the SMA CH3OH and 
13 CH 3 O H images (6U « 171 x 079, E up /k = 110- 730 K) ob- 
tained by Beut her et all (120051) around 337 GHz. It is interesting 
to note that dM-2 is close to the 25 GHz and 95 GHz CH 3OH 
maser emission peaks rep orted by iMatsakis et all (119801) and 
iPlambeck & Wrightl d!988l) . respectively. 



Population diagrams: We use population diagrams to estimate 
the column density of CH2DOH (e.g., [Goldsmith & Langer, 
fl999llTurnerl[T99Tl) 



In 



N up 3k f T B dV 
— " = ln ^ T~ 



Q 



kT ri 



(1) 



http://spec.jpl.nasa.gov/ 



where N is the total column density, T mt the rotational temper- 
ature, S the line str ength, and u the dipole m oment (jj. x = 1 .44 
D and fi z = 0.89 D. llvash & Dennisonl ll953h . In this work, the 
partition function Q is calculated according to our approxima- 
tion Q = 2.25 r 1 5 . We used the four CH 2 DOH lines at 223 
GHz (eo state) observed with the same spatial resolution in this 
population diagram analysis. Because our data points fall in a 
narrow energy range, the population diagrams give a tempera- 
ture of about 100 K with a very large uncertainty. Hence, the 
ro tational temperatu res were assumed and taken from the work 
of IFavre et al.1 (1201 ll) to reduce the statistical error. We assumed 
here that both CH3OH and HCOOCH3 trace the same gas be- 
cause both lines show the same emission peaks, Vlsr, and AV. 
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Fig. 1. Left: Channel maps of the CH2DOH 52,4 - 42.3 (E ap /k = 95.2 K) emission at 223422.3 MHz for a synthesized beam of 
l'/8 x(X'8. Contours run from 40 mJy beam -1 (2 cr) to 280 mJy beam -1 (4.84 K) in steps of 40 mJy beam -1 , and the dashed contours 
represent -20 mJy beam -1 . The bottom-right panel shows the integrated intensity (from 5 to 1 1 km s -1 ) in contours running from 
0.2 to 1.0 Jy beam -1 km s -1 in steps of 0.2 Jy beam -1 km s -1 , and the dashed contours represent -0.08 Jy beam -1 km s -1 . Middle: 
Channel maps of the CH3OD emission at 226185.9 MHz (E up /k = 37.3 K) for a synthesized beam of 3'.'6 x 2'.'3. Contours run from 
0.12 to 2.52 Jy beam -1 in steps of 0.04 Jy beam -1 (3 cr), and the dashed contours represent -0.08 Jy beam -1 . The bottom-right 
panel shows the integrated intensity of CH3OD (from 5 to 11 km s -1 ) in contours running from 10% to 90% in step of 10% of 
the peak intensity (3.1 Jy beam -1 km s -1 ), and the dashed contours represent -10% of the peak intensity. Right: Channel maps of 
the CH3OH 8_2 - 81 E (E up /k = 109.6 K) emission at 101469.8 MHz for a synthesized beam of 3'.'8 x 2'.'0. Contours run from 
40 to 490 mJy beam -1 (7.74 K) in steps of 50 mJy beam -1 . The dashed contours represent -30 (1 cr) and -90 mJy beam -1 . The 
bottom panel shows the integrated intensity of CH3OH (from 5 to 1 1 km s -1 ) in contours running from 0.3 to 3.0 Jy beam -1 km s -1 
in steps of 0.3 Jy beam -1 km s -1 , and t he d ashed contours represent -0.03 Jy beam -1 km s -1 . The black square marks the center 
of explosion according to lZapata et al.1 (120091) . The positions of source BN, source /, the Hot Core (HC), and IRc7 are marked by 
magenta triangles. The positions of deuterated methanol emission peaks (dM-1, dM-2, and dM-3) and KL-W are marked by red 
stars. Note that the channel separation in the CH3OD channel maps are resampled and only the three selected CH3OH channel maps 
are shown here for clarity (see Fig. IB.l I for the full CH3OH channel maps). 



Figure lB~3l shows the least-squares fit result where we assumed 
that the rotational temperatures according to the HCOOCH3 ex- 
citation temper atures are 130 K for dM-1 and 85 K for dM-2 
as derived bv lFavre et al.l(l201 ll) with the same spatial resolution 
(T/79 X 0'.'79). The derived CH 2 DOH column densities are given 
in Table 3. We find 8.8+0.9x 10 15 cm -2 and 2.4±0.3x 10 15 cm -2 
for dM-1 and dM-2, respectively. A dopting an H? col umn den- 
sity of 3 . 1 x 1 24 cm -2 as derived bv lFavre et al.l (1201 ll) from the 
223 GHz dust continuum, we obtain a CH2DOH relative abun- 
dance of 2.8 ± 0.3 x 10 -9 toward dM-1. The relative CH 2 DOH 
abundance is uncertain toward dM-2 because this clump does 



not have a strong dust continuum emission and no reliable H2 
column density can be derived. 

iJacq et all (1 19931) derived a CH 2 DOH column density of 
2.6 - 5.4 x 10 15 cm -2 toward Orion IRc2 with an estimated 
temperature of 88 K and an assumed source size of 15" from 
their IRAM 30m data. When smoothing our CH 2 DOH data to 
a source size of 15", the derived CH2DOH column density is 
0.2 - 1.4 x 10 15 cm -2 , assuming a maximum missing flux of 
50% and an excitation temperature of 88 K toward I Rc2. The 
column density disagrees with that of lJacq et al.l (1 19931) by a fac- 
tor of 2-3. Owing to the importance of deuteration ratios, which 
are obtained based mainly on previous single-dish estimates, we 
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Table 3. Derived column density and temperature for CH2DOH and CH3OH 



Position R.A. Dec. Size r rot (CH 2 DOH) jV(CH 2 DOH) r rol (CH 3 OH) AT(CH 3 OH) [CH 2 DOH]/[CH 3 OH] 

05 h 35 m -05°22' (K) (cm- 2 ) (K) (cm" 2 ) 



dM-1 


...14*442 


...34'.' 86 


V.'i 


X078 


130° 


8.8 + 0.9 x 10 15 
















J.'i 


x 2/0 


130" 


4.7 + 0.5 x 10 15 


130" 


4.2 


±0.4x 10 18 


1.1 ±0.2 x 10~ 3 








J.'i 


x 2'.'0 


60 


2.4 ±0.3 x 10 15 


60.3 ± 3.7 


2.1 


±0.3x 10 18 


1.1 ±0.2 x 10 3 


dM-2 


...14320 


...37'.'23 


V.'i 


xO.'8 


85" 


2.4 + 0.3 x 10 15 
















J.'i 


x 2'.'0 


85" 


1.7 + 0.1 x 10 15 


85" 


2.2 


±0.7x 10 18 


7.7 ± 2.5 x 10~ 4 








yn 


x 2'.'0 


66 


1.4 ± 0.2 x 10 15 


66.2 ± 4.7 


1.9 


±0.3x 10 18 


7.4 ± 1.6 x 10~ 4 


dM-3 


.. .14:107 


.. .37/43 


m 


xa/8 


58 


9.5 + 4.5 x 10 14 
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Notes. The [CH 2 DOH]/[CH 3 OH] ratios are derived only for the 3'.' 8 x 2'.'0 clump size. For dM-1 and dM-2, the CH 2 DOH column densities are 
derived from the population diagrams, and CH 3 OH column densities and temperatures are derived for all sources from the population diagrams. 
The CH 2 DOH column density of dM-3 is derived only from the 5 2/ t - 4 2 3 eo line at 223422.3 MHz, assuming a rotational temperature of 58 
K (estimated from CH 3 OH). The position s of dM-1, dM-2, dM-3, and KL-W are given at Epoch J2000.0. (a) The methyl formate rotational 
temperatures derived bv lFavre et al .1(1201 ll) are adopted here for CH 2 DOH and CH 3 OH. 



have revised the CH2DOH column density calculations together 
with CH3OD in Appendix[A]to compare them with our interfer- 
ometric r esults. We find th at the CH2DOH column density ob- 
tained bv lJacq et alj d 1993b is overestimated by a factor of 2-3, 
and the revisited value is thus consistent with our result men- 
tioned above. 

3.2. CH 3 OH 

Three E-type methanol lines around 101 GHz (see Table 2) have 
been detected; the line parameters toward different peaks are 
summarized in Table IB. 21 The CH3OH emission exhibits a V- 
shaped structure (Figs. fTland lBTt similar to that observed for 
several other molecular lines and for dust continuum emission. 
The eastern side of the V-shaped structure follows the dense 
ridge, and the bottom part coincides with the Compact Ridge 
region. The western part of this V-shaped structure lies in the N- 
S direction, an d one CH3OH emiss ion peak (Orion KL western 
clump, KL-W, [Wright et all [1991 MF-4/5, iFavre et al.L l201lh 
is located to the north close to the IR source IRc6, which lie s 
~ 3" north of IRc7 (see, e.g.. iGezaril [l992l iGezari et al.L Il998l) . 
Because of the lower spatial resolution (3'.'79 x l'/99) of these 
CH3OH lines, only dM-1, dM-3, and KL-W are spatially re- 
solved. The average FWHM line width toward the CH3OH emis- 
sion peaks is 4.4 + 1 .2 km s _1 , and the average line width at KL- 
W is about 6.5 km s , which is likely due to the blending of 
two or more velocit y componen t s at 8 and 10-11 km s _1 close 
to MF-4 and MF-5 dFavre et all 1201 lb . In addition, the E-type 
methanol lines' average velocity is 8.2 ±0.9 and 8.6 ±0.9 km s _I 
toward the three CH2DOH emission peaks and KL-W, respec- 
tively. Since KL-W is located ~ 7" from the Hot Core and ~ 10" 
from the Compact Ridge, the 8 km s velocity component ob- 
served in single-dish observations (6>mb > 20") may contain a 
large amount of gas (~ 25% of our CH3OH total flux) from the 
KL-W region or the western side of the V-shaped structure. This 
will affect the analysis based only on the velocity components. 

B ecause the CH3OH l ines are optically thick in most cases 
(e.g., iMenten et all 1 1988b . we first estimated the optical depth 
of the CH3OH lines detected here around 101 GHz. In the 
Rayleigh-Jeans approximation we derive the opacity from 



where 7\, and T ex are the source brightness and excitation tem- 
peratures (using the derived rotational temperatures of 40-70 K), 
7bg is the background emission temperature, including both the 
cosmic microwave background radiation (2.73 K) and the Orion 
dust emission (around 0.8-1.4 K from the PdBI data, correcting 
for a high interferometric filtering loss in continuum of 90%). 
The average CH3OH optical depth is estimated to be < 0.2 for 
the three 101 GHz CH3OH lines. Therefore, in the optically thin 
case, we can apply the population diagram method to estimate 
the E-type methanol column densities and rotational tempera- 
tures. The population diagrams are shown in Figure |B~4"1 and the 
derived excitation temperatures and column densities are listed 
in Table 

Rotational temperatures of 40-70 K and column densities 
of about 10 18 cirr 2 are found toward selected clumps in Orion 
BN/ KL, assuming a n A/E methanol abundance ratio of 1.2 
(IMenten et allll988b . The derived C H3OH column densiti es are 
higher than the IRAM 30m result of IMenten et all dl988b . who 
obtained a CH3OH column density of 3.4+0. 3x 10 17 cm" 2 with a 
rotational temperature of 128 + 10 K, including both narrow and 
broad components and assuming a source size of 25". To com- 
pare ours with their results, we smoothed our CH3OH spectra to 
the beam size (25") of the 30m telescope at 101-102 GHz (Fig. 
|B31 We derived a CH3OH column density of 2.3 ± 0.5 x 10 17 
cm -2 with a rot ation temperature of a bout 60 K at the same po- 
sition as used bv lMenten et alj d 1988b . Additionally, by adopting 
a higher temperature of 130 K, we derived a CH3OH column 
density of 4.6 + 0. 4 x 10 17 cm" 2 , whic h is close to the column 
density derived bv lMenten et all dl988b . 

Our methanol rotational temperature estimat e is lower than 
the 100-150 K obtained bv lMenten et alj d 19881) using th e 30m 
telescope and lower than the 80 K of iNeill et al.1 d201 lb using 
CARMA, but lies a t the low end of the 50-350 K obtained by 
Beuth er et al ] (120051) using SMA. Lower rotational temperatures 
may be due to the non-LTE behavior (or subthermal excitation, 
see § 14. 3t of methanol. In addition, choosing only the lower- £up 
transitions tends to result in a lo wer excitation temperature from 
popul ation diagrams (see, e.g., iBlake et ail 119871; iParise et all 
2002). Nevertheless, column densities are not very sensitive to 
rotational temperatures in this temperature range (Fig. IB. 61 , and 
the differences in CH3OH column densities derived with differ- 
ent temperatures are less than a factor of about 3. 
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3.3. CHiOD 

A new transition of CH3OD (5_i - 4_i E) was detected at 
226185.9 MHz for the first time in the ISM at the edge of our 
spectrometer and imaged toward Orion BN/KL. The CH3OD 
spectra for the selected sources in Orion BN/KL are shown in 
Figure [2] and the line parameters are listed in Table IB. 31 The 
CH3OD channel maps are shown in Figure [T] The CH3OD and 
CH2DOH spatial distributions are very similar (see Fig.[3]c), and 
both maps show the strongest emission at dM-1 with weaker 
emission toward dM-2 and dM-3. The line widths of CH3OD 
toward dM-1 and dM-2 are about 2-3 km s _1 , whereas two ve- 
locity components are clearly seen toward dM-3 with line widths 
of about 1-2 km s _1 . 

3.4. HDO 

In one of our low angular resolution data sets (3'.'6 x 2'.'3), we 
detected one HDO line at 225896.7 MHz (3 1>2 - 2 2 ,i). This is 
not sufficient for a detailed analysis, but it allows us to compare 
the spatial distribution and spectra between deuterated methanol 
and water emissions. Figure [2] clearly shows that the optically 
thin lines from rarer isotopologs of methanol have similar spec- 
tral profiles with an average line width of about 3 km s _1 and 
LSR velocities of about 8 km s _1 at dM-1, dM-2, and dM-3. On 
the other hand, the HDO spectra in Orion BN/KL have broader 
line widths, which are likely composed by two or more veloc- 
ity components. For example, the HDO spectra at dM-1 can be 
decomposed by two components at LSR velocities of about 5 
and 7 km s _! with line widths of about 10 and 4 km s _1 , respec- 
tively. The narrow line feature is likely produced by the same gas 
at which isotopologic methanol lines have been detected, and 
the broad line feature may be due mainly to the shock-heated 
gas associated with the Hot Core and source I region where the 
hot ammonia and highly excited CH3 CN are present (see, e.g., 
Godd i et all 1201 lHZapata et al.Ll20Tll) . 



4. Discussion 

4.1. The CH 2 DOH/CH 3 OH abundance ratios 

To investigate the CH3OH deuteration across the Orion BN/KL 
region, we smoothed our CH2DOH data to the spatial resolution 
of the E methanol lines (3'.'8 x 2'.'0) and derived the population 
diagram (Fig. IB.3b by using methanol rotational temperatures 
and HCOOCH3 rotational temperatures for low and high tem- 
perature reference, respectively (see Table [3] which summarizes 
our calculations toward different clumps in Orion BN-KL). The 
[CH 2 DOH]/[CH 3 OH] abundance ratios are 1.1 +0.2 x 10" 3 , and 
7.4 + 1.6 x 10~ 4 , and 8.8 + 0.8 x 10~ 4 for dM-1, dM-2, and dM-3 
with lower temperatures, respectively. An upper limit of about 
2 x 10~ 4 is given for the [CH 2 DOH]/[CH 3 OH] abundance ratio 
at KL-W. It is interesting to note that the [CH 2 DOH]/[CH 3 OH] 
ratio of dM-3 is higher than that of dM-2. 

These abundance ratios do not depend much on the as- 
sumed r iot (see Fig. IB. 61 ). If higher excitation temperatures 
are adopted (e.g., 130 K for dM-1 and 85 K for dM-2), the 
[CH2DOH]/[CH30H] ratios do not change much. These ra- 
tios (observed toward different clumps in Orion BN/KL, see 
Table [3]) are about one order of magnit u de low er than the ra- 
tio (0.01-0.08) reported by iJacq et al.1 (Il993l) . However, the 
[CH2DOH]/[CH30H] abundance ratios derived for the selected 
sources (0.8 - 1.3 x 10 3 ) are consistent with a revised ratio of 
0.8-2.8 xlO 3 , where possible line blending and different observ- 



ing positions were taken into account in the revisited calcula- 
tions of the 30m observations (see AppendixlAlfor more details). 
In addition, the [CH 2 DOH]/[CH 3 OH] ratios in Orion BN/KL 
are lower than the deuterated methanol and formaldehyde abun- 
dance ratios estimated in several low-mass protostars, e.g., up to 
about 0.6 for rCHz DOHMCHjOHl an d 0.3 for [HDCO]/[H 2 CO] 
in Class sources (Pai-is e et al.l l2006). 

Assuming that both CH 2 DOH and CH 3 OH are in LTE 
and have the same excitation temperature, we can produce 
a [CH 2 DOH]/[CH30H] ratio map by using two transitions 
of CH 2 DOH and CH3OH and taking into account their line 
strengths, partition functions, and upper- state en e rgies, the 
same method as used recently by Rataic zak et al.l d201 lb . In 
the [CH 2 DOH]/[CH 3 OH] ratio map (Fig. |3] d), dM-1 ex- 
hibits a higher ratio than observed toward the southern part 
of Orion BN/KL where dM-2 and dM-3 are located. It is in- 
teresting to note that the source / region also shows higher 
CH 2 DOH/CH30H intensity ratios on each side of the bipo- 
lar outflow (marked with a double black arrow in the image). 
Since the feature associated with source I is only seen in the 
strongest CH 2 DOH line at 223422.3 MHz, more observations 
are needed to confirm this result. Because we derive an upper 
limit of 2.2 X 10~ 4 for the [CH 2 DOH]/[CH 3 OH] abundance ratio 
toward KL-W, which is well below the ratios derived for dM-1, 
dM-2, and dM-3, there might be different physical and/or chem- 
ical conditions at KL-W. Nevertheless, we cannot exclude yet 
that the non-detection of the CH 2 DOH emission toward KL-W 
may result from spatial filtering in our PdBI interferometric data. 

The deuteration ratios derived here are not markedly differ- 
ent from one clump to another except perhaps for KL-W (where 
an upper limit is derived). However, it is clear that the deutera- 
tion ratio tends to decrease from the dense clumps (dM-1 to dM- 
3) to the more extended gas in Orion BN/KL. The H 2 column 
densities of dM-1 and dM-3 are about 3 - 5 x 10 24 cirT 2 , but 
the rotational temperature of HCOOCH3 toward dM-1 (about 
130 K) is higher than toward dM-2 and dM-3 (about 80 K). 
Therefore, the slight difference observed in the deuteration ratios 
is likely due to different gas temperatures of the clumps. Based 
on our present knowledge, we suggest that various competing 
processes could be at work. For instance, global, slow heating of 
the Orion BN/KL region through ongoing star formation should 
produce widespread effects and the expected net result is a fairly 
homogeneous deuterated ratio. On the other hand, heating by lu- 
minous infrared sources (e.g., in the vicinity of source 7) could 
result in a higher temperature in this region, leading to a higher 
deuteration ratio toward dM- 1 . 

4.2. The CH 2 DOH/CH 3 OD abundance ratios 

Although only one CH3OD line (5_[ - 4_i E at 226185.9 
MHz, E up /k = 37.3 K) has been detected in our observations, 
a similar [CH 2 DOH]/[CH30D] abundance ratio image can be 
produced in the same fashion as [CH 2 DOH]/[CH30H] by us- 
ing the CH 2 DOH 3 1,3 - 2 0j2 o line (E up /k = 35.7 K) at 
225878.3 MHz. Figure [3] (e) shows the [CH 2 DOH]/[CH 3 OD] 
abundance ratio image by assuming the same excitation tem- 
perature for CH 2 DOH and CH3OD, the method used for the 
[CH 2 DOH]/[CH30H] ratio map. This abundance ratio does not 
depend much on the temperature (see Fig. IB. 6 1 b). Our re- 
sults show that the [CH 2 DOH]/[CH 3 OD] abundance ratios are 
less than unity in the central region of Orion BN/KL and 
the mean ratio is 0.7 + 0.3, which also differs from the ra- 
tio (1.1-1.5) reported bv lJacq et al.l (119931) . However, the revis- 
ited calculation as mentioned above (Appendix lAl of the previ- 
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Fig. 2. Spectrum comparison of HDO, CH 3 OH, 13 CH 3 OH, CH 2 DOH, and CH 3 OD toward selected sources in Orion BN/KL with 
a similar resolution of about 3". The SMA 13 CH30H line (Y.-W. Tang, priv. comm.) is plotted in blue for a synthesized beam of 
3'.'0 x 2'.'0. The CH 3 OH line is plotted in red for a synthesized beam of 3'.'8 x 270. The HDO (thick black lines), CH3OD (thin black 
lines), and CH2DOH (magenta) lines are plotted for a synthesized beam of 3'.'6 x 2'.'3. Dotted lines indicate a Vlsr of 8 km s -1 . Note 
that the spectrum intensities of CH3OH and CH2DOH are multiplied by 2. Clearly, the spectra of methanol isotopologs have narrow 
line widths compared with those of HDO. 



ous 30m observations of C H 2 DOH dJaca et all 1 19931). CH3OD 
dMauersbergeret all 119881) . and CH3OH dMenten et all [l988) 
gives a lower [CH 2 DOH]/[CH 3 OD] abundance ratio of < 0.6 
toward IRc2 (adopting a source size of 15"), which is consis- 
tent with our results and is much lower t han the values derived 
in low -mass protostars (e.g. JParise et all 12002; Ratai czak et all 
120111) . Additionally, we note that the [CH 2 DOH]/[CH 3 OD] 
abundance ratios are higher at dM-1 and dM-3 compared with 
that at dM-2. 

The fact that the [CH 2 DOH]/[CH 3 OD] ratios are less than 
unity in the central Orion BN/KL region is surprising because 
most of the theoretical and experimental studies on deuterated 
meth anol predict more abundant CH?DOH than CH3OD (see , 



e.gJCharnlev et allll997tlOsamura et alll2004lNagaoka etail 
2005: lRataiczak et alll2009l) . Since CTLDOH and CH 3 OD have 
very similar emission distributions in Orion BN/KL, their for- 
matio ns seem to be closely related. As shown by lOsamura et al.l 
(12004 . exchanging protons and deuterons on the two dif- 
ferent parts of the methanol backbone is very inefficient in 
the post-evapo r ative gas-phase. However, the simulations of 
Osamura et al ] (12004 show that the water vapor that evapo- 
rates from the surface above 80 K strongly affects the abun- 
dance of CH3OD, because the protonation of H 2 and HDO in 
gas phase leads to additional fractionation throughout the reac- 



tion network. Therefore, the observed abundance ratio between 
CH 2 DOH and CH3OD might not reflect the primitive abundance 
ratio in the ice mantle of grain surfaces. This may explain the 
much lower [CH 2 DOH]/[CH 3 OD] ratios seen in Orion BN/KL, 
where most of the regions have higher temperatures than 80 K. 
Besides, it indicates that H 2 and HDO may be two critical 
molecules for understanding the formation processes of deuter- 
ated methanol in Orion BN/KL. 

4.3. CH^OH subthermal excitation? 



Recently, Wan g et all d201 ll) detected more than 300 methanol 
lines around 524 GHz (0 MB * 43") and 1061 GHz (6» MB ~ 20") 
with the HerschelfHlFl instrument. Although these data were 
obtained with a broad beamwidth, it may be useful to investi- 
gate the implications of the Wang et al. results on our own data 
in terms of CH3OH excitation. The narrow and broad CH3OH 
components observed by Wang et al. show average FWHM line 
widths of 2.7 and 4-11 km s -1 at V LS r = 8.6 and 6-9 km s -1 , 
respectively. Their LSR velocities are consistent with our data. 
However, their broad components' line widths are wider than 
ours. This is likely due to optical depth broadening (r ~ 15 in the 
Herschel/HIFI CH3OH data), although we cannot exclude that 
some faint broad and extended velocity components detected by 
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Fig. 3. (a) Color image showing the integrated intensity map of the CH3OH 8_2 - 81 E line (E up /k — 109.6 K, integrated from 5 
to 1 1 km s -1 ) for a synthesized beam of 3'.'8 x 2'.'0. Black contours represent the CH2DOH 52,4 - 42,3 line emission (E up /k = 48.4 
K, integrated from 5 to 11 km s~') and run from 15% to 95% in steps of 10% of the peak intensity (1.0 Jy beam -1 km s -1 ) for a 
synthesized beam of l'.'S x 078. Black dashed contours correspond to -5% of the peak intensity, (b) Color image showing the PdBI 
integrated intensity map of of the HDO 3i,2 - 22,1 line {E up /k = 167.6 K, integrated from -8 to +17 km s -1 ) for a synthesized beam 
of 3'.'6 x 273. Black contours represent the SMA 13 CH 3 OH 12 U i - 12 ,i2 A-+ emission (E up /k = 192.8 K, integrated from to 12 
km s _1 , Y.-W. Tang, priv. comm.), running from 20% to 90% in steps of 10% of the peak intensity (316.6 Jy beam -1 km s -1 ) for 
a synthesized beam of 370 x 270. Black dashed contours correspond to -10% of the peak intensity, (c) Integrated intensity map of 
the CH3OD 5_i - 4_i E line (E up /k = 37.3 K, from 5 to 11 km s -1 ) overlaid with the CH 2 DOH 3i, 2 - 2 , 2 o line (E up /k = 35.7 
K, integrated from 5 to 1 1 km s -1 ) in black contours for a synthesized beam of 376 x 273. Black contours run from 0.28 (2 <x) to 
1.12 Jy beam" 1 km s -1 in steps of 0.14 Jy beam -1 km s _1 , and black dashed contours correspond to -0.07 Jy beam -1 km s _1 . (d) 
[CH 2 DOH]/[CH 3 OH] abundance ratio map for a resolution of 378 x 270. The [CH 2 DOH]/[CH 3 OD] abundance ratio is calculated 
by assuming the same excitation temperature (120 K) for both molecules, (e) [CH2DOH]/[CH3(DD] abundance ratio map for a 
resolution of 376 x 273. The [CH2DOH]/[CH.30D] abundance ratio is calculated by assuming the same excitation temperature (120 
K) for both CH2DOH and CH^OD. Blue contours run from 0.4 to 1.0 in steps of 0.1. (f) The black square marks the center of 
explosion according to Zap ata et al.l (|2009). The positions of source BN, the Hot Core (HC), IRc6/7, and source //n, and SMA-1 
are marked by triangles. The positions of deuterated methanol emission peaks (dM-1, dM-2, and dM-3) and KL-W are marked by 
stars. The bipolar outflow of source / is indicated. 



Herschel may also be partly filtered out in our interferometric 
data. 

Our three CH3OH 101 GHz transitions include upper energy 
levels with E up /k < 1 10 K, from which we derived T rol =40-70 
K. The CH3OH lines in this e nergy level range detected by 
Herschel (see, e.g., Figure 4 in IWang et all l20lTl) also give a 
similar trend: a low rotational temperature and a high column 



density are derived by using only the low E ap data in the CH3OH 
population diagram. The CH3OH lines from the Herschel/HIFI 
observations are optically thick, but the 13 CH30H lines (t = 
0.03-0.28) also show the same trend. This suggests that our low- 
temperature result (T mt = 40-70 K) is unlikely to be caused by 
optical depth effects. Because the LSR velocities of our CH3OH 
lines are consistent with the Herschel/HIFI data, it is more likely 
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Fig. 4. Channel maps of the CH 2 DOH emission at 225878.3 MHz (E up /k = 35.7 K) in black contours overlaid with the HDO 
3i,2 - 22.i line (E up /k = 167.6 K) at 225896.7 MHz in light-blue contours for a synthesized beam of 3'.'6 x 273. The HDO emission 
contours run from 10% to 90% in steps of 10% of the peak intensity (17.9 Jy beam -1 ), and the dashed contours represent -0.2 
Jy beam 1 (5 cr). The CH2DOH emission contours run from 0.08 to 0.48 Jy beam 1 in steps of 0.08 Jy be am -1 , and the blac k 
dashed contours represent -0.04 Jy beam 1 (1 cr). The black square marks the explosion center according to Zap ata et al.l J2009). 
The positions of source BN, the Hot Core (HC), and IRc7 are marked by magenta triangles, and the positions of the deuterated 
methanol emission peaks (dM-1, dM-2, and dM-3) are marked by red stars. 



that our low CH3OH excitation temperatures are not caused by a 
cold component along the line of sight, but that they are largely 
generated by sub thermal excitation conditions. 

4.4. Comparison with the HDO emission 

In Figure |U the superposition of HDO 3 12 -22,1 and CH2DOH 
3 1,3 - 2q,2 lines shows a globally similar distribution for different 
velocity channels. This result is significant because both lines 
were detected in the same PDBI data set around 225.9 GHz. 
However, details differ, as one may expect for two transitions 
with different upper energy levels (168 K and 36 K for HDO and 
CH2DOH, respectively), and for different molecular formation 
scenarios, which imply, for instance, that methanol is formed on 
grain surfaces at higher densities than water. 

The HDO emission peaks betwe en the Orion Hot Cor e and 
dM-1, near the dust clump SMA-1 dBeufher et all 12004 ; see 
also our Fig. |3J. Dust emission toward this peak has been de - 
tected with SMA in the 870 fim band by iTang et all (l2010l) . 
Owing to a position offset of about 1", it is not clear whether 
this clump coincide s with the SMA-1 position observed by 
Beuth eretaT1(l2004l) with a resolution of 0778x0765. In addition 



to the prominent emission peak near SMA-1, the HDO emission 
shows a similar V-shaped distribution as CH3OH. There is HDO 
emission close to the IR source IRc7, but it is interesting to note 
that the HDO emission does not peak at IRc7. It peaks to the 
southwest of IRc7 instead. A similar situation is observed with 
SMA for the 13 CH30H emission (see the right panel of Figure 
|3j Y.-W. Tang, priv. comm.). In addition, it seems that the over- 
all CH2DOH emission distribution is shifted to the south by a 
few arcseconds with respect to the HDO emission. For exam- 
ple, the strongest HDO emission peak , which correspond s to the 
NH3 column density peak derived bv lGoddi et al.l (1201 ll) . is lo- 
cated 3" north of dM-1. This shift is clearly seen in Figure [3] 
where the HDO and I3 CH30H emissions lie farther north than 
the CH3OH and CH2DOH emissions. The different spatial dis- 
tributions observed for CH2DOH and CH3OH on one side, and 
HDO and 13 CH30H on the other side (see Figs. [3] and HJ cannot 
be attributed to optical depth effects because (a) we have shown 
that our CH2DOH (1.3 mm emission) and CH3OH (3 mm emis- 
sion) data are optically thin; (b) our estimate o f the HDO line 
opacity is < 0.3 (agrees wifh jjacq et all 1 1990b . The difference 
observed in the distribution of these molecular species could be 
due to different temperatures across the Orion molecular mate- 
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rial involving processes such as heating by luminous infrared 
sources or shocks. 



4.5. Deuteration ratios of water and methanol 

The [HDO]/[H20] abundance ratio has been reported to be 
about a few times 10~ 4 toward Orion BN/KL (see, e.g., 
iPardo et all 120011; iJacq et al.L 119901) . which is about one order 
of magnitude lower than the [CH 2 DOH]/[CH 3 OH] ratio derived 
here. [HDO]/[H 2 0] ratios lower than the [CH 2 DOH]/[CH 3 OH] 
ratios have also been reported t oward low-mass protostars 
dParise et all l2005t iLiu et all 120111) . Since heavy water (D 2 0) 
is much less abundant than HDO, i.e., rD 2 01/[HD01 = 1.7 x 10 3 
in IRAS 16293-2422 dButner et alll2007l) . HDO seems to be the 
main deuterium (gas phase) reservoir for water in Orion BN/KL. 
For methanol, CH3OD is as abundant as CH 2 DOH, and other 
doubly or multiply deuterated methanol species are expected 
to be less abundant. Therefore, the total deuteration ratio for 
methanol (including CH 2 DOH and CH3OD) in Orion BN/KL 
is about one order of magnitude higher than that of water (in- 
cluding HDO and D 2 0). If both water and methanol sublimate 
from ice mantles with similar time scales, and supposing that 
their destruction rates are comparable, the higher abundance of 
deuterated methanol may just reflect different chemical reac- 
tion speeds on dust surface. At a later stage, formation and de- 
struction rates in the gas phase will of course complicate the 
analysis. Some experiments have been conducted to investigate 
this question. For example, H-D exchange between water and 
methanol has been proposed by several authors in different phys- 
ical conditi ons, e.g., temperature increase and ultrav iolet light 
irradiation (iRataiczak et all 120091; IWeber et all 120091) . In addi- 
tion, the enrichment of deuterated methanol has been experimen- 
tally reproduced by H-D substitution in solid methanol at 10 K. 
However, a high atomic D/H ratio of 0. 1 was r equired and the 
substitution was not seen in water and ammonia (Nagaoka et al., 
120051) . Therefore, H-D substitution may qualitatively explain 
why methanol is easier to be deuterated than water; but more 
chemical experiments and modeling are clearly needed to quan- 
titatively address this question. 



5. Conclusions 

The main findings and conclusions of our study based on ob- 
servations of several transitions of deuterated methanol and one 
transition of deuterated water in Orion BN/KL are as follows. 

1 . We have obtained the first high angular resolution ( l'/8 x078) 
CH 2 DOH images detected around 223.5 GHz toward Orion 
BN/KL and compared these data with somewhat lower res- 
olution (3'.'8 x 270) CH3OH images at 101.5 GHz. The 
strongest CH 2 DOH and CH3OH emissions come from the 
Hot Core southwest region exhibiting an LSR velocity of 
about 8 km s _1 , typical of the Orion Compact Ridge re- 
gion. The CH 2 DOH emission is clumpy and the column 
densities are estimated to be about 1 - 9 x 10 15 cirT 2 
toward these clumps. The CH3OH column densities are 
about 3 - 5 x 10 17 cm" 2 across Orion BN/KL, leading to 
a [CH 2 DOH]/[CH 3 OH] deuteration ratio of 0.8 - 1.3 x 10 -3 
toward three deuterated methanol clumps and below 2 x 10 
toward KL-W. 

2. The [CH 2 DOH]/[CH 3 OD] abundance ratio map was ob- 
tained for Orion BN/KL, and their ratios are less than unity at 
the central part of the region. These ratios are lower than the 



statistical factor of 3 derived in the simplest deuteration mod- 
els, and definite ly lower than the values derived in low-m ass 
protostars (e.g. JParise et all 120021: IRataiczak et all 1201 11) . 

3. We have mapped with moderately high spatial resolution 
(376 x 273) the 225.9 GHz transition of HDO and compared 
its distribution with CH 2 DOH, CH 3 OH, and 13 CH 3 OH. We 
find that the deuterated water ratio is about one order of mag- 
nitude lower than the deuterated methanol ratio. H-D substi- 
tution may explain why methanol is easier to be deuterated 
than water. 

4. The deuteration ratios derived in this work are not strongly 
different from one clump to another, except perhaps toward 
KL-W where more observations are desirable to conclude. 
However, to explain the slight differences observed locally 
in the abundance ratios of identified clumps, we suggest that 
various processes could be competing, for instance, heating 
by luminous infrared sources, or heating by shocks. 
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Appendix A: Revisiting the CH 2 DOH/CH 3 OD 
abundance ratios in Orion BN/KL 

To compare our high angular resolution data with previous 
single-dish observations, we re-c alculated the temperatures 
and column densities obtained by Mauersberg er et al.l <[l988) 
and iJacq et al.l (1 19931) from their 30m data for CH 3 OD and 
CH2DOH, respectively. The line parameters of CH3OD and 
CH2DOH used in the calculations are listed in Table IA.1I and 
the population diagrams are shown in Figure IA.1I 

CH 3 OD 

We selected only five clean CH3OD lines (without any 
apparent line-blending) around 140-156 GHz detected by 
Mauersb erger et al.l (119881) with a similar resolution of 16 "—17" 
to reduce t he fi tting uncertainty. The CH3OD lines detected by 
UacqetalJd 19931) were excluded because no CH3OD spectra can 
be used to judge the possible line-blending and data quality. 
Effective line st r ength s S u 2 of the CH3OD lines were taken from 
lAnderson et al.l (119881) . and the partition function Q = 1.41 T 1 5 
was used. The rotational temperature (132.4±19.6 K) and col- 
umn density (4.4 + 0.3 x 10 15 cm -2 ) were estimated by a least- 
squares fit ("Fig. IA.lt . assuming a source size of 15". The revised 
CH3OD temperature and co l umn d ensity are consistent with 
those of iMauersberger et al.l (Il988l) . who derived a rotational 
temperature of 50-150 K and a column density of 1 - 5 x 10 15 

-2 

cm . 
CH 2 DOH 

Fo r CH2DOH we ju dge that only the three detections reported 
by IJacq et af] (1 19931) can be reliably used in the population di- 
agra ms (see Table | A. 11 1. Both old CH2DOH Sp 2 values taken 
from lJacq et al.l(ll993l) and new ones (see Table lA. lb provided by 
B. Parise (priv. comm.) were used in the calculations. (Fig. IA.ll i. 
The CH2DOH partition function Q = 2.25 T 1 5 and a source size 
of 15" were used here. 

The new values of CH2DOH Sp 2 em ployed in th i s work are 
based on those tabulated in the paper by P arise et al.l (120021) but 
a factor of three higher. This revision corre cts for an inconsis - 
tency with the partition function reported bv lParise et al.l(l2002l) . 
and some dipole moment values can be quite uncert ain, which 
may e xplain the discrepancy with the values used by [Jacq et al.l 
d 1993b and why they were never officially released in the JPL 
databas e. Howeve r , as sh own below, both the old Sp 2 values 
used bv lJacq et alJ (1 19931) and new ones adopted here give simi- 
lar estimates of the CH2DOH column densities in Orion BN/KL. 
We believe that the uncertainties of the CH2DOH Sp 2 values of 
our present work are not dominant in the column density cal- 
culations. Nevertheless, a new spectroscopic study of CH2DOH 
with a robust fitting of the Hamiltonian would definitely give 
better confidence in these intensities. 

Using the old Sp 2 , a least-squares fit to the CH2DOH data 
shows a low rotational temperature (about 25 + 10 K) with a col- 
umn de nsity of 5.0 X 10 14 cm 2 , consistent with the first esti- 
mate of IJacq et alJ ( 1 1993b . who derived a rotational temperature 
of 34^" K and a column density of 6x 1 14 cm' 2 , However, since 
the rotational temperatures derived by Menten et al. (1988J and 
IMauersberger et al.l (119881) fo r CH 3 OH and CH 3OD are about 
100 K toward Orion BN/KL. Ilacq et alJ (1 19931) also adopted a 
higher temperature for CH2DOH. They obtained a CH2DOH 
column density of 3.9 x 10 15 cm 2 , which is 2.6 times higher 
than our revised value (1.6 x 10 15 cm" 2 , see Fig. IA.1I upper 
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Fig.A.l. Population diagrams of CH3OD and CH 2 DOH us- 
ing th e 30m data of Mauersberge r et all (H988) and Lfacq et all 
( 1993 9, respecti vely. The old CH2DOH effective line strengths of 
Jacq et al. ( 1993) (upper panel) and the new ones (lower panel) 
were used in the calculations. The beam filling factor was taken 
into account by assuming a source size of 15". The rotational 
temperatures and column densities were estimated by a least- 
squares fit, and the results and uncertainties are given in the di- 
agrams. For CH2DOH, the red dashed line in the upper panel 
in dicates the fit with a fixed temperature of 100 K as adopted 
bv lJacq et al.l (Tl993) and the one in the lower panel the fit with 
a temperature of 130 K (the same as the CH3OD temperature 
derived here). 



panel). This is because they appear to correct the beam-filling 
factor twice in their Ni/gi calculations and populatio n diagram- 
Theref ore, the CH2DOH column density derived by IJacq et aLl 
(Il993l) with a temperature of 100 K was overestimated by a fac- 
tor of 2-3. 

With the new CH2DOH Sp 2 , we derive a somewhat higher 
rotational temperature (40.6 + 14.4 K) and a column density of 
5.0 x 10 14 cm" 2 similar to that derived by using the old Sp 2 . 
This suggests that the uncertainties in the CH2DOH effective line 
strengths do not strongly affect the CH2DOH column density 
estimate in this case. 

CH 2 DOH/CHt, OD abundance ratios 

The overestima ted CH2DOH column density derived by 
IJacq et all (1 19931) leads to a higher [CH 2 DOH]/[CH 3 OD] abun- 
dance ratio (1.1-1.5) toward IRc2. The new calculation indicates 
that the [CH 2 DOH]/[CH 3 OD] abundance ratio is 0.2-0.5 with 
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a source size of 15", assuming CH2DOH and CH3OD have the 
same temperature of about 130 K. For a lower CH2DOH temper- 
ature (i.e., about 40 K), the [CH 2 DOH]/[CH 3 OD] ratio is even 
lower (about 0. 1). 

In addition, the CH3OD data of Mauersbe rger et all (1 19881) 
were taken a t the p osition close to KL-W (KL-Wm hereafter, 
Men ten etall l"l988) instead of IRc2, which was used as the cen- 
tral position by iJacq et all (119931) . According to our CH3OD 
image (Fig. |3j, the intensity observed toward IRc2 is about 
twice as high as that of KL-Wm, assuming a 30m telescope 
beam size of about 18" at 140 GHz and the same filtering. 
Th erefore, the CH3OD abun dance in Orion BN/KL derived 
bv iMauersber geret al.l (1 1988 ) may be underestimated by 50% 
at most. A similar calculation (Fig. \B.5\ also shows that the 
CH3OH column density toward KL-Wm is about 20% lower than 
that toward IRc2. 

In short, we conclude that the revisited 
[CH2DOH]/[CH30D] abundance ratio derived from the previ- 
ous 30m observations is < 0.6, and the [CH 2 DOH]/[CH 3 OH] 
abundance ratio is estimated to be 0.8 - 2.8 x 10~ 3 toward 
IRc2, by adopting a CH3O H column density of 4.7 + 0.3 x 10 17 
cm -2 dMenten et all Il988l) and correcting the possible density 
underestimate of about 20% toward IRc2 (for a source size of 
15"). 



Appendix B: Complementary figures and tables 
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Table A.l. 30m CH 3 OD and CH 2 DOH line parameters 



Frequency 


Transition 


Eup/k 




J T MB dV 


r 




Nulgtf Comment 
(cm" 2 ) 


(MHz) 


(Ska 


,0 


(K) 


(D 2 ) 


(Kkms- 1 ) 






CH3OD " 


140175.20 


4l,3 


- 4 ,4 A— 


19.8 


9.3 


5.1 ±0.7 


0.42 


1.55 ±0.25 x 10 12 


143741.65 


5 1 .4 


- 5 ,5 A — 


27.5 


11.2 


6.6 ±0.7 


0.43 


1.57 


± 0.20 x 10 12 


148359.77 


60,6 


- 5 15 A++ 


31.7 


5.7 


3.3 ±0.3 


0.45 


1.45 


±0.16x 10 12 


153324.00 


7l,6 


- 7 ,7 A— 


47.4 


14.7 


7.6 ± 0.6 


0.47 


1.21 


±0.11 x 10 12 


155533.08 


ll.o 


- ,o E 


7.8 


1.0 


0.9 ± 0.4 


0.47 


1.93 


± 1.09 x 10 12 


CH 2 DOH" 


99672.23 


61.5 


- 60,6 e 


50.4 


12.3" 


0.3 ±0.1 


0.27 


1.51 


±0.61 x 10" 


136151.26 


3 1,2 


- 2i,i e a 


17.1 


2.1" 


0.6 ±0.1 


0.41 


8.58 


±1.72x10" * 


226818.36 


5l,4 


- 4i,3 e 


36.7 


3.8" 


1.6 ±0.2 


0.66 


4.72 


±0.71x10" * 


99672.23 


61,5 


- 60,6 e 


50.4 


5.8' 


0.3 ±0.1 


0.27 


3.23 


± 1.29 x 10" 


136151.26 


3 1,2 


- 2i,i e 


17.1 


3.1' 


0.6 ±0.1 


0.41 


5.91 


±1.18x10" * 


226818.36 


5 1,4 


- 4|, 3 e 


36.7 


5.5'' 


1.6 ±0.2 


0.66 


3.27 


±0.49x10" 



Notes. <fl) iMauersbereer et all (1 19881) : centered at KL-W M (05 b 35 m 14 s .18, -05°22'26'.'5, J2000). {b) Uacq et alj il993l) : centered at IRc2 
(05 h 35 m 14 s .47, -05°22'3a'2, J2000). (c) The beam filling factor / = <9 s 2 /(<9 2 + 6> 2 ) is calculated by assuming a source size 8 S of 15". 6 b is the 
30m beam size. (< " The beam filling factor is taken into account in th e calculation of N u , and the rms noise and calibration uncertainties are in- 
cluded (assuming 20%). {e) The Old Sfi 2 taken from lJacq et al.1 i ll 9931) ( " The New S/j 2 provided by B. Parise (priv. comm.). <s) The data obtained 
at the offset position (0",-6") of IRc2. 



Table B.l. PdBI CH 2 DOH line parameters at dM-1, dM-2, and dM-3 



Frequency" 


Transition 




sn 2 




Vlsr 


AV 


J TdV 


Comment 


(MHz) 




(K) 


(D 2 ) 


(K) 


(km s" 1 ) 


(km s- 1 ) 


(Kkms- 1 ) 




dM-1 


105806.4100 


Hi, 10 - llo.n e\ 


159.4 


6.11 


0.45 ± 0.02 


7.81 ±0.89 


2.85 ± 1.09 


1.23 ±0.05 


blended c 


223422.2629 


52,4 - 4 2 ,3 e 


48.4 


5.15 


4.77 ± 0.35 


7.96 ± 0.84 


2.71 ± 1.14 


14.29 ± 0.50 




223616.1420 


54,2 _ 4 4 l , 54,1 - 4 4 , eo 


95.2 


2.15* 


2.91 ±0.35 


7.70 ± 0.84 


2.78 ±2.14 


6.95 ± 0.23 




223691.5380 


53,3 - 4 3 ,2 e 


68.2 


3.96 


4.81 ±0.35 


6.99 ± 0.84 


2.78 ± 1.68 


15.87 ± 1.76 


blended'' 


223697.1880 


5i,2 - 4 3 ,i e 


68.2 


3.96 


3.04 ± 0.35 


7.79 ± 0.84 


1.64 ± 1.35 


6.17 ± 1.03 


blended c 


225878.2540 


3 1,3 _ 2o,2 Co 


35.7 


2.28 


1.27 ±0.12 


7.75 ± 0.42 


2.67 ± 0.44 


3.10 ±0.20 




dM-2 


105806.4100 


1 1 1,10 _ 1 lo.ii e \ 


159.4 


6.11 


0.38 ± 0.02 


7.81 ±0.89 


2.91 ± 1.07 


1.02 ±0.05 




223422.2629 


5 2 ,4 - 42,3 eo 


48.4 


5.15 


2.79 ± 0.35 


7.96 ± 0.84 


1.83 ±0.54 


5.19 ±0.48 




223616.1420 


54,2 - 44,1 , 54,1 - 44,o «0 


95.2 


2.15* 


1.22 ±0.35 


7.30 ± 0.84 


2.84 ± 1.33 


3.03 ± 0.32 




223691.5380 


53,3 - 4 3 ,2 eo 


68.2 


3.96 


1.40 ±0.35 


7.79 ± 0.84 


2.80 ± 1.68 


2.43 ± 0.55 


blended rf 


223697.1880 


53,2 -43,i e 


68.2 


3.96 


0.68 ± 0.35 


7.86 ± 0.84 


1.22 ±0.97 


1.89 ±0.55 




225878.2540 


3 1,3 - 2q,2 Oq 


35.7 


2.28 


1.07 ±0.12 


7.76 ± 0.42 


2.52 ± 0.42 


2.29 ±0.19 




dM-3 


105806.4100 


1 1 1,10 - 1 lo.ii e\ 


159.4 


6.11 


0.32 ± 0.02 


7.89 ± 0.89 


3.36 ±0.89 


1.02 ± 0.05 




223422.2629 


52,4 - 4 2 , 3 e 


48.4 


5.15 


0.76 ± 0.35 


7.92 ± 0.84 




2.87 ± 1.35 




223616.1420 


54,2 _ 4 4 ,i , 54,1 - 4 4 ,o e 


95.2 


2.15* 


< 0.35 






< 1.75 




223691.5380 


53,3 - 4 3 ,2 e 


68.2 


3.96 


1.40 ±0.35 


7.76 ± 0.84 


3.43 ± 0.90 


4.23 ± 0.26 




223697.1880 


53,2 - 4 3J e 


68.2 


3.96 


0.93 ± 0.35 


7.84 ± 0.84 


4.18 ± 1.05 


3.43 ±0.31 




225878.2540 


3 1,3 - 2o,2 00 


35.7 


2.28 


0.76 ±0.12 


7.71 ±0.42 


0.87 ± 0.45 


0.71 ±0.11 


1st component 










0.66 ±0.12 


9.05 ± 0.42 


1.22 ±0.45 


0.83 ±0.18 


2nd component 



Notes. <D) Line frequencies of the detected lines {b) Effective line strengths of each single line {c) This line is blended with the H 13 CCCN J = 12-11 
line at 105799.1 130 MHz toward dM-1. {d) This line is blended with the U223694.8 line. (c) This line is blended with the U223694.3 line and the 
acetone line at 223692.1 MHz. 
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Fig.B.l. Left: Channel maps of the CH 2 DOH doublet (E up /k = 48.4 K) emission at 223616.1 MHz for a synthesized beam of 
l'/8 x(X'8. Contours run from 40 mJy beam -1 (2 cr) to 160 mJy beam -1 (2.77 K) in steps of 40 mJy beam -1 , and the dashed contours 
represent -20 mJy beam -1 . The bottom-right panel shows the integrated intensity (from 5 to 1 1 km s -1 ) in contours running from 
0.15 to 0.45 Jy beam -1 km s -1 in steps of 0.15 Jy beam" 1 km s , and the dashed contours represent -0.08 Jy beam -1 km s -1 . Right: 
Channel maps of the CH 3 OH 7_ 2 -7i E (E np /k = 91 .0 K) emission at 101293.4 MHz for a synthesized beam of 378 xZ'O. Contours 
run from 50 to 350 mJy beam -1 in steps of 50 mJy beam -1 . The dashed contours represent -30 mJy beam -1 (l cr). The bottom-right 
panel shows the integrated intensity of CH3OH (from 5 to 1 1 km s -1 ) in contours running from 0.3 to 1.5 Jy beam -1 km s -1 in steps 
of 0.3 Jy bea m -1 km s" 1 , and the dashed contours represent -0.03 Jy beam -1 km s -1 . The black square marks the center of explosion 
according to ZaDat aetaLl (120091) . The positions of source BN, the Hot Core (HC), IRc7, and source / are marked by triangles. The 
positions of deuterated methanol emission peaks (dM-1, dM-2, and dM-3) and KL-W are marked by red stars. 
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Fig. B.2. CH 2 DOH spectra toward dM-1, dM-2, and dM-3 for a synthesized beam of 1'.'8 x (X'8. The dashed lines indicate a V LS r of 
8 km s _1 . The unidentified molecular lines are denoted with their frequencies, (a) For dM-1, the CH2DOH 1 1 110 - 1 ln,n e\ line is 
blended with a broad H 13 CCCN J = 12-11 line (AV ~ lOkms 1 ) at 105799.1 MHz. Therms noise (1 &) is 0.02 K.'(b)-(d) The 
dotted line indicates the CH 2 DOH 5 3>3 - 4 3>2 e Q line at 223691.5 MHz, which is blended with the (CH 3 ) 2 CO 17 7j n - 16 6 ,io EA/AE 
lines at 223692.1 MHz. The rms noise (1 a) is 0.35 K. (e) The strong line to the left is the HDO 3i, 2 - 2 2 ,i line at 225896.7 MHz. 
Two velocity components at 7.8 and 9.0 km s _1 are seen toward dM-3. The rms noise (1 <x) is 0.12 K. 
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Fig. B.3. CH 2 DOH population diagrams of dM-1 and dM-2. (a)-(b) Population diagrams of the high angular resolution data. The 
rotational temperatures are assumed to be 130 K and 85 K for dM-1 and dM-2, respectively, (c)-(d) Population diagrams for the 
data smoothed to a resolution of 3'.'8 x 2'.'0. Low rotational temperatures and higher ones are assumed to be 60 K and 130 K for 
dM-1 and 66 K and 85 K for dM-2, respectively. The fits of low temperatures are plotted in black lines and higher temperatures in 
blue dashed lines. The CH 2 DOH column densities are estimated by a least-squares fit, where the uncertainties include the statistical 
error, the rms noise, and calibration uncertainties (10%). 
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Fig. B.4. CH3OH spectra (3'.'8 x 2'.'0) toward the four positions identified in Fig. [T]in the Orion BN/KL region. The methanol 
rotational temperatures and column densities are estimated by a least-squares fit, where the uncertainties include the statistical error, 
the rms noise, and calibration uncertainties (10%). Blue dashed lines indicate the fits with fixed temperatures. Black dashed lines 
indicate a Vlsr of 8 km s _1 . 
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Fig. B.5. Simila r plot as Fig.IB.4l CH tOH spectra were smoothed to the 30m beam size of 25" to compare with the previous 30m 
observations by M enten et al. ( 19881) . Upper panels show t he CH ^O H spectra and popula t ion di agram taken at the same position 
(05 h 35 m 14 s .18, -05°22'26'.'5, J2000) used bv lMenten etal] (fl988h and lMauersbereeret al.1 JT988l) . which is located about 5" to the 
north of KL-W (denoted as K L-Wm)- Lower pa nels show the CH3OH spectra and population diagram taken at IRc2 (05 h 35 m 14 ; !47, 
-05°22'30'.'2, J2000) used bv lJacq et alJ dl993h . Blue dashed lines indicate the fits with a fixed temperature of 130 K. Black dashed 
lines indicate a Vlsr of 8 km s~'. Note that the CH3OH column density derived toward KL-Wm is about 20% lower than that of 
IRc2. 
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Fig. B.6. (a) Variation of derived CH3OH column densities with respect to different rotational temperatures toward dM- 1 . The corre- 
sponding [CH2DOH]/[CH30H] abundance (derived in population diagrams) ratios are shown in red. The best-fit result corresponds 
to a CH3OH column density of 2.1 x 10 18 cirT 2 and a temperature of about 60 K. For temperatures lower than 200 K, the differences 
in the derived column densities are clearly smaller than a factor of about 3. The [CH2DOH]/[CH30H] abundance ratios derived 
from different temperatures are lower than 1 .1 x 10~ 3 . (b) Variation of the relative abundance (based on the ratio of two transitions 
of CH2DOH and CH3OD as seen in Fig.[3]e) of the two isotopologs toward dM-1 as a function of the adopted temperatures. Note 
that this ratio does not depend much on the temperature. 



Table B.2. PdBI CH3OH line parameters at dM-1, dM-2, dM-3, and KL-W 



Frequency 


Transition 


E up /k 


V 


* peak 


Vlsr 


AV 


J TdV 


(MHz) 




(Jk) 




(K) 


(D 2 ) 


(K) 


(km 


r 1 ) 


(km 


r 1 ) 


(Kkms- 1 ) 












dM 


-1 (05 h 35 ra 14!442, -05° 


22' 34'.' 86) 










101185.4530 


6_ 


-2 


-6, 


E 


74.7 


0.021 


3.72 ± 0.05 


7.01 ± 


0.93 


4.09 ± 


1.58 


17.92 


±5.56 


101293.4150 


7. 


-2 


"7! 


E 


91.0 


0.046 


6.50 ± 0.05 


8.31 ± 


0.93 


4.29 ± 


0.95 


30.72 


±5.61 


101469.8050 


8_ 


-2 


-81 


E 


109.6 


0.091 


8.59 ± 0.05 


7.84 ± 


0.93 


4.18 ± 


0.65 


43.49 


±5.45 












dM 


-2 (05 h 35 ra 14!320, -05° 


22' 37'.' 23) 










101185.4530 


6_ 


-2 


-6, 


E 


74.7 


0.021 


2.70 ± 0.05 


7.86 ± 


0.93 


4.33 ± 


1.98 


14.48 


±5.38 


101293.4150 


7. 


-2 


"7! 


E 


91.0 


0.046 


5.29 ± 0.05 


8.26 ± 


0.93 


4.48 ± 


1.19 


25.49 


±5.52 


101469.8050 


8_ 


-2 


-81 


E 


109.6 


0.091 


7.36 ± 0.05 


8.76 ± 


0.93 


4.40 ± 


0.75 


37.07 


±5.34 












dM 


-3 (05 h 35 m 14!107, -05° 


22' 37'.' 43) 










101185.4530 


6_ 


-2 


-6, 


E 


74.7 


0.021 


2.98 ± 0.05 


8.81 ± 


0.93 


4.46 ± 


1.82 


14.28 


±4.65 


101293.4150 


7. 


-2 


"7! 


E 


91.0 


0.046 


5.00 ± 0.05 


8.31 ± 


0.93 


4.59 ± 


1.14 


23.78 


±4.76 


101469.8050 


8_ 


-2 


-81 


E 


109.6 


0.091 


6.84 ± 0.05 


8.76 ± 


0.93 


4.54 ± 


0.74 


33.54 


±4.62 












KL- 


W (05 h 35 


n 14 s .159, -05 


22'28'.'25) 










101185.4530 


6. 


-2 


-61 


E 


74.7 


0.021 


2.70 ± 0.05 


8.76 ± 


0.93 


6.52 ± 


6.13 


8.65 ± 6.05 


101293.4150 


7_ 


-2 


"7! 


E 


91.0 


0.046 


4.16 ±0.05 


8.31 ± 


0.93 


6.72 ± 


3.44 


14.09 


±5.65 


101469.8050 


8. 


-2 


"8! 


E 


109.6 


0.091 


5.40 ± 0.05 


8.76 ± 


0.93 


6.17 ± 


2.24 


17.38 


±5.25 



Notes. The R. A. and Dec. coordinates at the J2000.0 epoch are given for the four sources. Frequency, E uv jk, and 5 fi 2 data were taken from the JPL 
database. The Vlsr is measured at the peak temperatures, and the line widths and integrated intensities are estimated by fitting a Gaussian profile. 



Table B.3. PdBI CH3OD 5_i - 4_i E line parameters at dM-1, dM-2, and dM-3 



Position 


^peak 


Vlsr 


AV 


J TdV 


Comment 




(K) 


(km s- 1 ) 


(km s" 1 ) 


(Kkms"') 




dM-1 


1.91 ±0.24 


7.69 ± 0.42 


2.19 ±0.42 


4.47 ± 0.27 




dM-2 


1.82 ±0.22 


7.58 ± 0.42 


2.87 ± 0.42 


5.55 ±0.36 




dM-3 


1.60 ±0.30 


7.43 ± 0.42 


1.39 ±0.42 


2.36 ±0.21 


1st component 




0.83 ±0.27 


9.28 ± 0.42 


1.94 ±0.41 


1.71 ±0.28 


2nd component 



